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Abstract 



We generalise Gaussian/Wuppertal smearing in order to produce non-spherical wave functions. 
We show that we can achieve a reduction in the noise-to-signal ratio for correlation functions of 
certain hadrons at non-zero momentum, while at the same time preserving a good projection 
on the ground state. 



1 Introduction 



In lattice QCD, masses, energies and transition matrix elements of asymptotic states with 
given quantum numbers are typically extracted from the Euclidean time dependence of suit- 
able correlation functions. The contribution of the lowest state can be separated from those 
of other states by inserting the source fields at large-enough time distances. The associated 
statistical error can be estimated from the spectral properties of the theory [TJ|2] and usually 
grows exponentially with the time separation^ In practice, this makes it very difficult to find 
a window where statistical and systematic errors are both under control. The use of smeared 
sources (i.e. interpolating fields with small overlaps on the excited states) ameliorates the sit- 
uation, as the contribution from the lowest state is then expected to saturate the correlation 
function already at small Euclidean separations. 

In general, the problem becomes more severe as non-zero spatial momenta are considered, 
as for instance in the computation of form-factors, and it is therefore worth trying to optimise 
the smearing procedure depending on the kinematics. The most prominent example is given 
by the pion. At zero momentum the noise-to-signal ratio (-Rns) is expected to be constant as 
a function of the time xq, whilst for p / 0, from a generalization of the arguments presented 
in [TJ[2] one expects 

R m (x ) OC e W™l+V?-rn«)x ^ for X ^ OO . (1) 

Throughout the paper we will characterise the momentum of the pion (nucleon) by a vector 
of integers, q = (91,(72,93), where p = qp In figure [TJ we show the effective energies 
extracted from two-point correlation functions of fields interpolating pions at non-zero spatial 
momentum q for different values of q 2 . The plots refer to actual simulation results from 
measurements performed on an ensemble of L/a = 32, T/a = 64, Nf = 2 unquenched gauge 
configurations generated within the CLS initiative using 2 flavours of non-perturbatively 0(a) 
improved Wilson fermions. Further numerical details can be found in [3j[4]. In total, 168 
configurations were used and for each of the configurations four evenly spaced source positions 
were utilised to improve the available statistics. We implemented Wuppertal smearing 5,6 
at the source and sink with parameters k = 2.9 and n = 140, whose meaning will be explained 
in the next section. These values were chosen to maximise the length of the effective mass 
plateau (i.e. p 2 = 0) for pions and nucleons and produce a Gaussian wavefunction with a 
radius of approximately 0.5 fm |2j. 

The dramatic deterioration of the signal as q 2 increases is quite striking. In figure § the 
function log i?Ns( 2; o) from the same dataset is plotted as a function of xq for different values of 
q 2 and it is compared to the expected asymptotic behaviour given in eq. ([!]). The asymptotic 
trend sets in quite early in the example and it is consistent with the noise being dominated 
by a zero-momentum two pion state. 

In the following, we will show how a significant improvement can be obtained by adopting 
non-spherical quark-smearing in order to interpolate boosted hadrons, instead of using the 
spherically symmetric function optimised for the zero-momentum it is usually done. 

One exception is the pion channel (at zero momentum) , for which the signal-to-noise ratio is asymptotically 
constant as a function of the source-sink separation. 
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Figure 1: Effective energies for pions at different values of q 2 

We underline that whilst this method provides a reduction of the noise in the measured 
quantities (the effective energies in this study) , it does not achieve an exponential improvement 
of the noise-to-signal ratio. 

A preliminary account of our work has been given in [7j . Here we study both the pion and 
nucleon channels. Prom the comparison it is clear that a much bigger gain in noise reduction 
for boosted pions can be achieved due to the fact that, for non-zero lattice momenta, the 
Lorentz factor 7 is much larger for pions than it is for nucleons. 



2 Anisotropic smearing for mesonic and baryonic 2-point functions 

We start by recalling the construction of Gaussian/Wuppertal smearing |5,6|. Given the 
hopping operator in direction k 

H k (x,x>) = U k {x)5 x , x+% + Ul(x - k)8 x , x _ k , (2) 

the smeared fermion field ip^ is defined as 

^)(x) = (l + KH) n iP(x) , (3) 

with 

H = H X + H y + H z . (4) 
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Figure 2: Noise-to-signal ratio for different q plotted against the asymptotic behaviour of eq. il 



In our generalization, we simply promote K and H to spatial- vectors k and H with components 
x, y and z and define 

^(«)(x) = (1 + K-H) n ^(x) . (5) 

In this study, we fix n = 140 and smear the gauge links in the hopping operator by apply- 
ing one level of hypercubic smearing [8j. For the particular choice of smearing parameters, we 
follow the definition of "HYP2" links in |9j. We will consider two-point correlation functions 
of the type 

C(x) = (o(^\x),^ K \x))O\^\0)^ {K \0))) , (6) 

where the operator O interpolates mesons or baryons and is therefore bi- or tri-linear in the 
fields ip^ and ip^. Our goal is to reduce the zero- momentum component, which domi- 
nates the noise, or equivalently to enlarge the high-momentum components in the four-point 
correlation function describing the variance. This is achieved by spatially squeezing the quark- 
wavefunction in the momentum direction. The choice is intuitively motivated by the physical 
picture that boosted hadrons are squeezed in the momentum direction because of the rela- 
tivistic space-contraction. As long as the size of the resulting wavefunction in the different 
directions is comparable to that of a boosted hadron we do not expect to spoil the good 
overlap with the ground state and therefore the fast approach to a plateau in the effective 
energy. 

This generalization of Gaussian/Wuppertal smearing allows us to mix and parametrically 
deform point-like wavefunctions and broad Gaussian ones by tuning the components of k in 
eq. ([5]). The first is known to give a better signal, especially at non-zero momentum, whereas 
the latter results in a better plateau 10 . We show in figure [3j the yz-section of a spherical 
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wavefunction corresponding to 2.9 (left panel) as well as the ?/2:-section of an 

anisotropic one, obtained using k x = n y = 2.9, k z = 2.9/16 (right panel). The lattice used is 
32 3 x 64 points with a lattice spacing of 0.063(2) fm and a pion mass of roughly 450 MeV 3|ll . 
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Figure 3: The yz-section of a spherical wavefunction 2.9 (left panel) and an anisotropic 

wavefunction k x = n y = 2.9, k z = 2.9/16 (right panel). 

The anisotropic smearing is used so that the spatial squeezing is applied only to the 
directions with momenta. A numerical determination in one dimension demonstrated that 
for fixed n in eq. the size of the wave function in direction I depends on m as roughly 
k 1 / 5 . Figure [4] compares the effect of spherical and anisotropic smearing on the effective 
energies for the first four Fourier momenta of the pion. This demonstrates a clear gain in 
the statistical signal, whilst preserving a good overlap with the ground state. The gain is 
more pronounced for an increasingly boosted pion. For instance, for q = (0, 1, 1) and (0, 0, 2) 
the identification of the plateau and the statistical error is significantly improved. The plots 
display the direction of squeezing as well as the corresponding values of 7 to indicate the 
amount of boosting, where 



for a hadron H. 

The amount by which we spatially squeeze the directions with momenta should be de- 
pendent on the amount of boosting, 7. We checked this behaviour through implementing 
several choices for the wavefunction scaling 1/N, with N = 2, 4, 8, 16 and 32. Figure [4] 
demonstrates that the improvement in statistical noise has saturated by N = 16 and that any 
further improvement for TV = 32 is accompanied by a slightly reduced suppression of excited 
states. In order to further corroborate our intuitive argument we also checked that squeezing 
the z-direction for q = (1,0,0) or (1,1,0), i.e. squeezing in an orthogonal direction to the 
momentum, does not provide a gain. 

The gain achieved by anisotropic smearing is further illustrated in figure [5j which shows 
the energies extracted from plateau fits to figure |4j The fit-ranges are chosen to be the same 
for spherical and anisotropic smearing, and are determined through a x 2 /dof minimisation 
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Figure 4: Effective energies for the pion for different values of q 2 using spherical (blue) and anisotropic 
(red and green) smearing. 



with a maximal plateau length. The continuum dispersion relation, determined from the 
physical pion mass, is overlaid to indicate the expected behaviour. 

We apply the same treatment used for the pion to the nucleon, figures [6] & [7j However, 
the signal-to-noise ratio is worse to begin with and so we are restricted from looking at 
significantly boosted nucleons. This means 7 ~ 1 and therefore the gain is not as striking as 
for the pion. The most likely explanation for the observed difference between the dispersion 
relation and the data points is the presence of cutoff effects. The latter can be estimated 
from a comparison between the continuum dispersion relation and a linear fit of the three 
data points with the smallest momenta for both the spherical and anisotropic smearing. The 
difference of 4% amounts to less than two sigma and is therefore not statistically significant. 
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3 Conclusions and outlook 

We have provided numerical evidence that by a simple generalization of Gaussian/Wuppertal 
smearing, which allows us to produce non-spherical quark-wavefuctions, the noise-to-signal 
ratio in correlation functions describing the propagation of hadrons at non-zero momentum 
can be significantly reduced. The prescription yields a reduction of the statistical uncertainty 
without affecting the onset of a plateau in the effective energy and consequently the overlap of 
the wavefunction on the ground state. The strongest evidence is for the case of pions, as those 
are already highly boosted at small values of the typical lattice momenta. The comparison 
between pions and nucleons also indicates that the statistical gain mainly depends on the 
boosting factor, 7, rather than on the mass of the hadron. For this reason we believe this 
method to be even more effective as the chiral limit is approached. Of course, this expectation 
requires further studies, which we will consider in the future. 

Recently, a very similar approach has been proposed in fl2| for nucleon interpolating 
fields. In that case, only a small improvement has been observed, which is consistent with our 
findings. 

The largest gain can probably be achieved by combining the approach with variational 
techniques, as in the Generalized Eigenvalue Problem method [13[[l4| . As a different smearing 



should in principle be adopted for each momentum, the use of all-to-all [15, 16 stochastic 
propagators could be advantageous in order to contain the numerical cost by reducing the 
number of inversions of the Dirac operator to be performed. We will explore these possibilities 
in future studies. 
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Figure 6: Effective energies for the nucleon for different values of q 2 using spherical (blue) and 
anisotropic (red and green) smearing. 

The main applications concern the determination of mesonic and baryonic form factors, 
e.g. for the B — > irlv process or for the nucleon form factors Ge and Gm, on which we have 



recently presented preliminary results 17 
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Figure 7: The dispersion relation plotted against the extracted effective energies for spherical (blue) 
and anisotropic (red) smearing for the nucleon. 



[2 
[3 

[4 

[5 

[6 

[s; 

[9 
[10 



G. P. Lepage, TASI 89 Summer School, Boulder, CO, Jun 4-30, 1989. 

P. Fritzsch, F. Knechtli, B. Leder, M. Marinkovic, S. Schaefer, R. Sommer and F. Vi- 
rotta, Nucl. Phys. B 865 (2012) 397 j |arXiv: 1205.53801 [hep-lat]]. 

S. Capitani, M. Delia Morte, G. von Hippel, B. Jager, A. Jiittner, B. Knippschild, 

H. B. Meyer and H. Wittig, |arXiv:1205.0"l80] [hep-lat]. 

S. Giisken, U. Low, K. H. Mutter, R. Sommer, A. Patel and K. Schilling, Phys. Lett. B 
227 (1989) 266. 

S. Giisken, Nucl. Phys. Proc. Suppl. 17 (1990) 361. 

M. Delia Morte, B. Jager, T. Rae and H. Wittig, "Improved interpolating fields for 
hadrons at non-zero momentum", PoS LATTICE 2012 (2012) 260, poster presented by 
BJ and TR at "Lattice 2012", Cairns, Australia, June 26, 2012. 



A. Hasenfratz and F. Knechtli, Phys. Rev. D 64 (2001) 034504 |hep-lat/0103029|. 



M. Delia Morte, A. Shindler and R. Sommer, JHEP 0508 (2005) 051 | hep-lat /0506008] 



H.-W. Lin, S. D. Cohen, R. G. Edwards, K. Orginos and D. G. Richards, arXiv:1005.0799 
[hep-lat] . 



[11] S. Capitani, M. Delia Morte, G. von Hippel, B. Knippschild and H. Wittig, PoS LAT- 
TICE 2011 (2011) 145 [arXiv:1110.6365l [hep-lat]]. 



8 



[12] D. S. Roberts, W. Kamleh, D. B. Leinweber, M. S. Mahbub and B. J. Menadue, 



arXiv: 1206.5891 [hep-lat]. 



[13] M. Liischer and U. Wolff, Nucl. Phys. B 339 (1990) 222. 

[14] B. Blossier, M. Delia Morte, G. von Hippel, T. Mendes and R. Sommer, JHEP 0904 
(2009) 094 |arXiv:0902.1265| [hep-lat]]. 

[15] J. Foley, K. Jimmy Juge, A. O'Cais, M. Peardon, S. M. Ryan and J. -I. Skullerud, 
Comput. Phys. Commun. 172 (2005) 145 [hep-lat/0505023] . 



[16] E. Endress, A. Jiittner and H. Wittig, |arXiv: 1111 .59 88 [hep-lat]. 

[17] S. Capitani, M. Delia Morte, G. von Hippel, B. Jager, B. Knippschild, H. B. Meyer 
and H. Wittig, "Excited state systematics in extracting nucleon electromagnetic form 
factors," PoS LATTICE 2012 (2012) 177, talk presented by TR at "Lattice 2012", 
Cairns, Australia, June 28, 2012. 



9 



